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Calculation of Propagation Parameters 
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An algorithm for computing the electromagnetic propagation modes 
and their associated propagation constants (i.e., the loss and phase 
shift per length) is rigorously developed for uniform cables. The con- 
ductors (including the shield, if present) are assumed to have circular 
cross sections and to be covered by two layers of dielectric insulation. 
By means of the algorithm, it is now possible to compute the propaga- 
tion parameters of uniform cables for all frequencies below the micro- 
wave range. The algorithm consists of calculating the eigenvectors and 
eigenvalues of a certain matrix, which are the set of conductor voltages 
of the modes and the associated propagation constants, respectively. 
The matrix in question is computed by methods developed in a com- 
panion paper on charge densities. 

I. INTRODUCTION 

Multipair cables consist of a collection of insulated wires surrounded 
by other dielectric materials, all of which are usually enclosed in a 
jacketed, metallic shield for mechanical protection and electrical isola- 
tion. The wires in a cable are generally helically twisted in pairs and, in 
practice, other, often undesired, nonuniformities occur along the cable. 
Nevertheless, the uniform cable, whose wires are straight and parallel, 
defined by the property that all its longitudinal cross sections are iden- 
tical, has successfully modeled some aspects of electrical propagation 
(such as loss per length) over a pair in a cable. 1 Also, some cables are 
directly modeled as uniform cables; so there is a need for studying uni- 
form cables. Moreover, further experience using the uniform-cable model 
may lead to a more exact model for nonuniform multipair cables. 

A rigorous analysis of electromagnetic propagation over uniform ca- 
bles, starting from Maxwell's equations, was performed by Carson 2 for 
cables having homogeneous dielectric material separating the wires. 
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Subsequently, Mead 3 applied his method to the shielded balanced pair 
to get a set of formulas for the primary constants (R, L,C,G). 

Later, Kuznetsov 4 developed another procedure for calculating the 
propagation parameters of uniform cables. He too assumed the dielectric 
separating the wires was homogeneous; also, he did not consider the 
presence of a shield. His procedure was abstractly stated, but he did 
derive a concrete formula for the loss and phase-shift per length in the 
special case of a single pair in free space (see Ref. 4, Chap. 2). Both 
Kuznetsov and Carson assumed that the frequency of excitation was high 
enough that the associated skin depth was less than the radius of the 
wires. 

In this paper, the electrical behavior of uniform cables is analyzed from 
first principles along the lines developed by Kuznetsov. The analysis 
pertains to all frequencies of excitation below the microwave range, it 
pertains to cables with or without a shield, and it provides for two layers 
of dielectric on the conductors* (including the shield). The final result, 
making use of results from the accompanying paper, 5 is an algorithm, 
suitable for use on a digital computer, that computes the various prop- 
agation parameters of the cable. 

The analysis starts with the concept of a modal solution. For a sinus- 
oidal excitation of frequency a>/27r, a modal solution for the electric field 
has the form, 

E = exp (iut - yz) E, 

where the Cartesian coordinate z runs along the cable, y denotes a 
propagation constant of the cable (giving the loss and phase-shift per 
length of the mode), and the vector field E is independent of time (t) and 
z. Thus, the electric field distribution for a modal solution is the same 
for all transverse cross sections of the cable, apart from a multiplying 
factor that gives the loss and phase shift per length for the mode, and 
similarly for the magnetic field. The general solution is a linear combi- 
nation of modal solutions. 

For propagating solutions, in particular, Kuznetsov has shown that 
there is a potential function V(i.e., a solution of Laplace's equation) 
defined in the dielectric-region, independent of z, and constant on the 
surface of the conductors, such that external to the conductors 

E= -VV + jX 

where i z denotes the unit vector in the 2 -direction and the function E z 
denotes the 2 -component of E. (His argument for this electrostatics 
approximation is enlarged in the next section to account for the shield 



* This permits the modeling of a new type of wire insulation recently introduced into 
cable manufacture: a layer of expanded insulation surrounded by a "skin" of solid insu- 
lation. The insulation on the shield is included for formal completeness. 
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and inhomogeneities in the dielectric). The various propagation modes 
of the cable, called here the excitation modes, are specified by the set 
of values for V on the conductors, which are simply the conductor volt- 
ages. 

The set of wire voltages at z = (the shield is always assumed to be 
grounded) is represented by the vector 

v = {vi, • • • , vm), 

where M is the number of wires in the cable. Certain particular volt- 
age-vectors correspond to excitation modes of the cable. Usually, there 
are M such normalized voltage vectors, each corresponding to a distinct 
propagation constant. But it can happen that two or more linearly in- 
dependent voltage vectors correspond to the same propagation constant; 
in this case, the propagation constant corresponds to the subspace 
spanned by these vectors, and the dimension of the subspace is called 
the multiplicity of the propagation constant. For either situation the 
collection of excitation modes and their associated propagation constants 
are the specific propagation parameters that are sought. The usual pri- 
mary constants of a transmission line or mode, if desired, can be deter- 
mined by standard formulas 6 from the propagation constant in con- 
junction with the capacitance matrix. 

Under the idealized conditions that the conductors are perfectly 
conducting and the dielectric is homogeneous, there is a single propa- 
gation constant, 

7 = ik e = ^M 1 ' 2 , (1) 

where 6 and n denote the permittivity and permeability of the dielectric, 
respectively. In such a case, the multiplicity of ik e is M. This is the so- 
called TEM mode, 6 where the longitudinal components of the electric 
and magnetic fields are zero. Other modes are evanescent, since the 
frequency range considered here lies below their cutoff frequencies. Since 
the dielectric in the cable may be inhomogeneous and the conductors 
are good but not perfect, the propagation solutions or modes are per- 
turbations of the TEM mode, the result being in general a set of excitation 
modes having distinct propagation constants. 

In the next section, the algorithm for computing the excitation modes 
and associated propagation constants is derived. The modes correspond 
to the eigenvectors of an M X M matrix, the propagation constants are 
simple functions of the associated eigenvalues, and the matrix itself is 
related to the capacitance and admittance matrices of the cable. In the 
third section, the algorithm is applied to the shielded balanced pair and 
a high-frequency approximation is indicated. Experimental testing of 
these results is reported in Ref. 7. 
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II. PROPAGATION MODES FOR UNIFORM CABLES 

In this section, an algorithm for computing the propagation modes 
and their associated propagation constants is derived for a uniform cable. 
The analysis pertains to cables having M straight and parallel wires of 
circular cross section, enclosed by a circular metallic shield. 

The geometry, materials, and excitation of the structure to be con- 
sidered are specified as follows: 

(i) The radius of the /nth conductor is denoted R m for m = 0, 1, •••, 
M, where fl refers to the inside radius of the shield. The shield has a 
uniform thickness th. 

(ii) The conductors have two layers of insulation with thickness (R m i 
- R m ) and (fi m2 -fi m i) for the wires (m ■ 1, -, M) and (R - R 0i ) and 
(Roi - Rq 2 ) for the shield. The permittivities (possibly complex) are e m 
and t m \ for the first and second layers, respectively, for m - 0, 1, -, M. 
In general, the subscript "0" will refer to the shield. 

(i«) The materials are nonmagnetic, and the permeability is n 
throughout. 

(iu) The material in the interstitial space outside the conductors and 
their insulation has permittivity e (possibly complex). 

(v) The excitation frequency co/2tt satisfies (a)/27r)/?n(«M) 1/2 « 1, that 
is, the wavelength is much greater than the radius of the cable. Since /? () 
generally exceeds 0.1 inch, this means that the frequency is below the 
microwave range. 

(ui) The conductivity of a conductor is denoted o m and we « a m for 

m = 0, 1, -, M. 

(vii) When the cross-sectional plane of the cable is viewed as the 
complex plane, the centers of the conductors are associated with the 
complex numbers b = Q,b\, —, 6m- A typical cross section is indicated 
in Fig. 1. 

The electrical behavior of a cable is described by the electric and 
magnetic fields, E and H, which satisfy Maxwell's equations. A modal 
solution, as discussed in Section I, has the form, 

E = exp(ia>£ — yz )E 

H = exp(iiot - yz)H, (2) 

where the vector fields E and H are independent of 2 and t. In terms of 
their longitudinal and transverse parts, 

E = E tr + i z E 2 

H = H tr + i z H z , (3) 

where the vector fields E lr and H lr are the projections of E and H, re- 
spectively, onto the transverse plane. Maxwell's equations imply that 

600 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1977 




Fig. 1 — Typical cable cross section. 



(see Ref. 4, page 3) 

-X 2 E ,r = yVE z + ia>ti VH Z X i z 
-X 2 H tr = yVH z + (k 2 /i<on) VE Z X i z , 
where taking ? as the permittivity as a function of position, 

in rath conductor 



k 2 = 



I k 2 m = -iwn a m m - 0, 1, • • • , M 
k 2 , = o) 2 ne outside the conductors, 



(4) 



(5) 



where the subscript "e" refers to the region exterior to the conductors, 



and 



X 2 = k 2 + y 2 . 



(6) 



The form of the solution in the space separating the conductors is 
determined by the electrostatics approximation. Kuznetsov 4 has vali- 
dated this approximation in three steps: 

(i) Both E zv and H Z( , and their derivatives are determined to better 
than one part in 2A/k P Ro as solutions of Laplace's equation. 

(«) There is a function V such that E' r = -VV and V 2 V = 0. 

(Hi) On each conductor the tangential part of E tr is zero; so V is 
constant on the conductors. 

For point (i), the original proof must be enlarged to deal with the pres- 
ence of a shield and inhomogeneities in the dielectric; this is done in 
Appendix A. Point («) is a direct consequence of point (0 because in eq. 
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(4) both the divergence and curl of VH Z X i z are zero, which implies it 
is the gradient of a solution to Laplace's equation, and point (Hi) comes 
from the estimate, {k 2 Jk 2 m « 1 for all m) [see Ref. 4, eq. (2.3) and Ap- 
pendix 1]. 

The potential which is one on the fcth conductor and zero on the others 
is denoted Vk k = 1, — , M. These form a basis set since for any potential 
V (which is zero on the shield) there are coefficients Uk k = 1,—, M such 
that 

M 

V=Y,VkV k . (7) 

Likewise, when E ze (k) denotes the longitudinal component that cor- 
responds to the potential V^ k = 1, —, M, then 

M 

E 2e = L u k E ze (k) (8) 

corresponds to the potential V. The coefficients Vk k = 1, — , M form the 
excitation vector 

v = (vi, ••-, vm)- (9) 

To find those excitation vectors that correspond to propagation modes 
and their associated propagation constants, a set of constraints is de- 
veloped, one for each conductor, by relating V and E ze at the conductor 
surfaces. At the outside surface of the mth wire ( W m ), eq. (4) implies 

„, „ dEze iun dHze , ,, /1A * 

-XlE Pme = y — + —— m=l,...M, (10) 

where (p m , <t>m) denote polar coordinates based at the center of the mth 
conductor. ButE Pme = -dV/dp m , so 

xi C"—^ = y Cir" d *™ m-i,-,M, (id 

JO dp m JO Opm 

where the derivatives are evaluated on W m . By eqs. (7) and (8), this is 
expressed as 

Xi £ »k \ ~ — d<p m = y £ v k I — d<p m 

k=l JO Op m k=\ JO Op m 

m=l,---,M. (12) 
The mk- element of the M X M capacitance matrix C is 



J»2tt $y k 
- d<j> m , 
OPm 



dp, 
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(13) 



and by eq. (4) for the transverse magnetic field, the mk- element of the 
M X M admittance matrix Y is 

r«-f w *,. »>*-**?■ r^^d^. (u) 

*S W m X- JO dp m 

Thus, the set of constraints in eq. (12) is jointly represented by the matrix 
equation, 

iu>Cv = yYv, (15) 

where the excitation vector v is viewed here as a column vector. 
But, as will be shown, the matrix Y has the form, 

y=|??, net 

X~ c 

where the matrix Y is independent of 7. Consequently, eq. (15) can be 
expressed 

Cv = (y 2 /X 2 e ) Yv or C" 1 Yv = (X 2 /y 2 ) v. (17) 

This means that to correspond to an excitation mode, the excitation 
vector v must be an eigenvector of the matrix C~ l Y, and conversely. If 
A is the associated eigenvalue, then the corresponding propagation 
constant 7 is 

7 = *M1 " A)~ 1/2 . (18) 

The remaining problem is to calculate the matrices C and Y. 

The elements of the capacitance matrix C can be computed by the 
method developed in Ref. 5 (see in particular Section III). This consists 
of inverting a matrix and is a standard operation on a digital computer. 
The elements of the Y- matrix are evaluated by deriving a boundary- 
value problem for the function E ze (in the region separating the con- 
ductors), then again applying methods from Ref. 5 to solve it. 

As the first step in determining Y, in Appendix B it is shown that E ze 
satisfies the following conditions: 

(i) V 2 E ze = 0, where V 2 denotes the Laplacian, (d 2 /dx-) + 

(dVdV 2 ). 

(ii) Across dielectric interfaces E ze is continuous and 5(dE ze /dp) = 
y8(dV/dp), where d/dp denotes the normal derivative and 8 denotes the 
difference or jump. 

(IB) (dE ze /d Pm ) - Dj: ze = y(dV/dp m ) at p m = R m for m = 0, 1, ~, 
M, where D w is a linear operator such that for m ~ 0,1, •••, M 

Du . e in<t, m _ (X cm /a nm ) e in <"n n=0, ±1, ±2, . . . , (19) 

and the coefficients a nm (which depend on a;) are defined in eq. (58) for 
m = and in eq. (52) for m 5* 0, both in Appendix B. In other words, in 
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condition (Hi), D^ alters the Fourier series of E ze on each conductor by 
multiplying the respective Fourier components. 
This boundary problem is simplified by putting 

E ze = yV + E' ze . (20) 

Then the corresponding conditions on the function E ze are 

(i) V*E' Mf = 0. 

(ii) Both E' ze and its derivatives are continuous across dielectric in- 
terfaces. 

(Hi) [bE'.Jbpm) - DJZ' ze = yD w V at p m = R m for m = 0, 1, •«, M. 

The function E' ZP restricted to a conductor surface is represented by 
a Fourier series expansion and this by an infinite vector whose compo- 
nents are the Fourier coefficients. If, consistent with the notation in Ref. 
5, u denotes the concatenation of these vectors for each conductor and 
p the concatenation of the vectors corresponding to e m (dE ze /dp m ) m 
= 0, 1, -,M, then condition (Hi) can be expressed by the matrix equa- 
tion, 

p - D^u = D a V, (21) 

where D^ denotes the diagonal matrix with elements t m X m /a nm for 
^m ^M and n — 0, 1, 2, •••. Also, there is a vector $ and matrices To, Ho, 
and G (defined in Ref. 5) such that 

u = T$ and p = G /3 - H u, (22) 

where the subscript "0" emphasizes the context of homogeneous di- 
electric, in accordance with condition (ii) for E' ze . It follows that 

(Go - H T - D a To) fl = t^V (23) 



or 



= -y[(H + DJT - Go]" 1 D W V. (24) 

But it follows from eq. (14) and eq. (16) that for 1 ^ m, k ^ M, 

- € m R m r2" dE' ze (k) ., , odd/ o a I 

Y m k ~ C mk = I — d<t> m = 2TrR m P ok /y = 2irefiok/y, 

T J ° dp - (25) 

where the latter equality is shown in Ref. 5, eq. (9). Therefore, the mk 
element of the f-matrix is computed by carrying out the matrix inver- 
sion in eq. (24) for V = V m to give (3 k- Of course, in practice, this in- 
version is performed on a truncated approximation of the matrix [(Ho 
+ DJTo-G ]. 

III. APPLICATION AND EXTENSIONS 

In this section, the general results of the last section are applied to the 
so-called shielded balanced pair. Also, extensions of the results to cables 
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without a shield and cables with circular holes in the interstitial dielectric 
are discussed. 

The shielded balanced pair (SBP) consists of a pair of identical wires 
placed symmetrically inside a circular shield (the center of the shield 
lies on the line joining the centers of the wires and midway between 
them). Since the structure is symmetric about the perpendicular bisector 
of the line joining the wire centers, there will be a balanced mode: the 
potential and fields will be antisymmetric about this line and, in par- 
ticular, the wire-voltages will be equal and opposite in sign. The problem 
is to determine the propagation constant (7) of the SBP. For convenience, 
the dielectric will be assumed to be homogeneous with permittivity e. 

The excitation vector v = (1, —1) is an eigenvector for the matrix C~ l Y 
and the corresponding eigenvalue is 

\=(Yu-Y l2 )/(C n -Cn). (26) 

Equation (18) gives 7 in terms of X. When V12 = V\ — Vo, 

dV V2 
'W\ dpi 

which is the total charge on wire no. 1 ( W\), and 



C n -C 12 =e(f ^ds = Q,, (27) 

Crwx dpi 



Yn-Y v2 =(e/y) (f ^ds. (28) 

These two quantities are evaluated as in section II (relative to the po- 
tential function VY2): Q\ by methods from Ref. 5, and Y n — Y v > by eqs. 
(24) and (25). 

A suggestive formula for 7 is obtained by applying Green's identity 
to eq. (28). Since both V12 and E ze are antisymmetric functions, and since 
V\2 — on the surface of the shield (Sh), 

<f ^*-i (f^Vnds, (29) 

%r W] dpi 2 kT v drj 

where Y denotes the entire boundary of the dielectric (W\, W±, Sh) and 
d/dr) denotes the normal derivative on T (into the dielectric). Since both 
E ze and V\-> satisfy Laplace's equation, Green's identity (see Ref. 8, Vol. 
2, page 252) implies that the derivative in the integrand can be switched 
to give 

+ - <f qoE ze ds, (30) 
2 tJrsh 
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where q\ and qo denote the charge densities on W\ and Sh, respectively. 
Thus, 

y = ik e \ l-(Qi)- 1 (f qi (E ze /y)ds 

-(2QJ-1 J? sh q (E ze /y)ds\ ; (31) 

the Fourier components of (E ze /y) — V12 on W, and Sh are To/3 with j8 
given by eq. (24). This formula for 7 indicates the contribution from each 
conductor. A comparison between measurement and theory for the SBP 
based on eq. (31) is presented in Ref. 7. 

For high frequencies (i.e., w for which conaRi » 1), the formula for E ze 
can be simplified by neglecting the first term of condition (Hi) (see Ref. 
4, eq 5.8, where X e should be X,). This gives 

E ze =-(y/e)D: l q m (m = 0, 1, 2). (32) 

This approximation together with the approximation 

(1 - 8)~ l/2 sl + %5 
gives 

7 - ik e fl + (2Q1)- 1 \$ 9iU>Z l Qi)ds 

V-2 fl sh Qo(D: ] qo)ds~\\. (33) 



+ 



The integrals in (31) and (33) can be expressed in terms of the Fourier 
coefficients of the charge density for convenience in calculating these 
expressions. 

The development in Section II can be extended to cables without a 
shield and to cables with circular holes in the interstitial dielectric. The 
calculation of the charge components, and hence, the capacitance and 
admittance matrices, in both of these situations was treated in Ref. 5, 
Section III. Again, eq. (17) must be solved for eigenvectors and eigen- 
values; but in the absence of a shield, the capacitance matrix is not in- 
vertible and the problem is what Kato 9 calls an eigen-problem in the 
generalized sense. In principle, these can be solved on a digital computer 
as before. 

The results presented here in conjunction with the computation 
techniques of the companion paper on charge densities allow compu- 
tation of the excitation modes and their associated propagation constants 
for uniform cables under a broad range of conditions. These include 
capacitive and resistive unbalance, caged shields that consist of a col- 
lection of unconnected wires, and many decades of excitation frequen- 
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cies. As the simplest, nontrivial example, the formulation was applied 
to the shielded balanced pair. As indicated in Ref. 7, the calculations 
agree reasonably well with measurements. 



APPENDIX A 

The Electrostatics Approximation 

The proof of the electrostatics approximation is broken into three 
parts in Section III. In this appendix the first part is verified; for the other 
two parts the reader is referred to Kuznetsov's work. 4 It is shown in 
particular that within an error for E ze of one part in (2A/k e Ro) (typically 
250 for 10 MHz), 

V 2 E 2e = 0, (34) 

and the same argument will hold for H ze . 
As indicated in eq. (1.5) of Ref. 5, 

(V 2 + X;,)E zl , = 0. (35) 

If E[\, is a function such that E° 2C = E ze on I\ E% is continuous across 
dielectric interfaces and b {dE^/dp) = 5(dE ze /dp) there, and 

V>E° ze = 0, (36) 

then E ze = E z \. + E\ v , where El,, = on T, El,, and its derivatives are 
continuous across dielectric interfaces, and 

(V 2 + Xl)E\„ = -X*E { l v . (37) 

When \px denotes the first eigenfunction of -V 2 that is zero on I\ then 
the L 2 -norm of E\,\s bounded by 

XI 



El. 



EU, (38) 



Xp Af 

where Xi is the eigenvalue associated with \pi. The smallest eigenvalue 
of —V 2 in the circle of radius Rq relative to a boundary condition of zero 
on the circumference is 2.4048//?o [2.4048 is the first zero of J (x)] and 
(see Ref. 8, Vol. 1, page 409, Theorem 3) 

Xy ^ (2.4048/flo) 2 - (39) 

Hence, 

\\E l J\ ^ (X e ny2.4048) 2 ||£° ( ,||, (40) 

and since X e ~ k e , 

E ze ^ E° ze , (41) 
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which was to be proved. Clearly, the same argument holds for the x and 
y derivatives of E ze . 

APPENDIX B 

Boundary Problem for E ze 

From the various conditions on E z , a boundary problem for E ze can 
be deduced. The result is stated as three conditions in Section II. In this 
appendix these conditions are verified. 

(i) From the electrostatics approximation, V 2 E ze = 0. 

(«) Since E z is a tangential component of E, it is continuous across 
interfaces. The difference in normal-derivative across the dielectric in- 
terfaces is deduced from eq. (4). Since H ze is continuous across inter- 
faces, 

-8(X^E pme ) = y8(^) m = 0,l,..-,M (42) 

at the interfaces. But by the usual boundary condition for the normal 
E field, 5(k 2 e E„ e ) = 0; so for m ■ 0, 1, •••, M, 

at the dielectric interfaces. 

(Hi) By the same reasoning as in (ii) 

dE ze dE zm 



dpm &Pr 



= -y(E pe - E pm ) atpm = R m , (44) 



where E zm and E pm refer to components of E in the mth conductor for 
m — 0, 1, —, M. Since 

E fjm =(k 2 Jk 2 m )E pe «E pe (45) 

for all m, eq. (44) to a good approximation becomes 

dE ze dE zm dV 

=7- — form = 0,l,--.,m. (46) 

<>Pm Op m dp m 

Finally, it is proved that 
dE 



'zm 



= D^mE^ for m = 0, 1, ■ • • , M, (47) 

?>Pm 

and since E zm = E ze at p m = R m , this will verify the third condition. This 
is proved first for the wires, then for the shield. 
Inside a wire, 

V 2 E zm + XlE zm = m=l,..-,M, (48) 

608 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1977 



and since y 2 « X; l} , the quantity X~ m can be replaced by k „, The solution 
for E zm has the form, 

fl«m"Sfl» 77*^ «rtfo»). (49) 

where the a„ are the Fourier coefficients of E z at p m = R m and J n denotes 
the rcth order Bessel function of the first kind. The radial derivative at 

Pm - Rm is 

— = k m L a /" ^ m * m) exv(in<fi m ) = k m D um E 2m , (50) 

where D^ is a linear operator, which for m = 1,—, M, 
D u , m --exp{in<t> m )^(k m /ot nm )exp(in<f> rn ) n = 0, ±1, ±2, • • • , (51) 
and 

«/n(&m#/7i) 



4* n \RmKm ) 



(52) 



Inside the shield, 

V* E z0 + HE z0 = 0, (53) 

where X,i has been replaced by k'i The solution has the form, 

E„ = Ya b ' ,J " ikm)) + CnYAkoP<)) expdn^), (54) 

where the a„ are the Fourier coefficients of E z at p () = R lh the 6„ and c„ 
are determined by reference to the boundary condition on the outside 
surface of the shield (p = Ro + th= R' ), and Y„ denotes the rcth order 
Bessel function of the second kind. Assuming the cable is not in the vi- 
cinity of sources or other conductors, the potential outside will be zero; 
hence, there is no charge on the outer surface of the shield, and by eq. 
(4) dE z /dpo is continuous at p = Ro> Since outside the cable, E z has the 
form 



Eg - £ fn (—) exp(w^o), 



/^ " exp(m0„), (55) 

for some /„, it follows from the continuity of dE z /dp () that in eq. (54), 



b n = Y n -i(koR' ) and c„ = -J n -i(k Ro) n - 0, ±1, ±2, • • • . (56) 
The radial derivative of E z at po = i?o is 

"7^ = feoE (aJctno) exp(m0„), (57) 

opo « 
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where 

_ Y n -i(koR' )J n (koRo) - Jn-i(koR' )Y n (koR ) 

an ° Yn-dkMJnikoRo) ~ Jn-l(koR'o)Y'n (kM 

n = 0, ±1, ±2, • • •. 
Hence, at po = Ro 

^ = koD uQ E z o, (59) 

where D w q is a linear operator as in eq. (51) with m = 0. The direct sum 10 
of the Du m for m = 0, 1, •«, M is D^, as given in eq. (19). 
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